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Abstract 

We consider the cosmological constraints on the holographic dark energy model 
by using the data set available from the type la supernovae (SNIa), CMB 
and BAO observations. The constrained parameters are critical to deter- 
mine the quintessence or quintom character the model. The SNIa and joint 
SNIa+CMB+BAO analysis give the best-fit results for j3 with priors on f2 m o 
and loq. Using montecarlo we obtained the best-fit values for f3, f2 m o and u>$. 
The statefinder and Ora diagnosis have been used to characterize the model 
over other DE models. 

PACS: 98.80. -k, 95.36,+x 

1 Introduction 

The astrophysical data from distant la supernovae observations pQ, [2], cosmic mi- 
crowave background anisotropy [3j, and large scale galaxy surveys pi], [5], all indicate 
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that the current Universe is not only expanding, it is accelerating due to some kind of 
negative-pressure form of matter known as dark energy ([E],[7J,[E],[H]). The combined 
analysis of cosmological observations also suggests that the universe is spatially flat, 
and consists of about ~ 1/3 of dark matter (the known baryonic and nonbaryonic 
dark matter), distributed in clustered structures (galaxies, clusters of galaxies, etc.) 
and ~ 2/3 of homogeneously distributed (unclustered) dark energy with negative 
pressure. Despite the high percentage of the dark energy component, its nature as 
well as its cosmological origin remain unknown at present and a wide variety of mod- 
els have been proposed to explain the nature of the dark energy and the accelerated 
expansion (see [HI 13 El IS] for review). Among the different models of dark energy, the 
holographic dark energy approach is quite interesting as it incorporates some concepts 
of the quantum gravity known as the holographic principle ( JTUJ, [TTJ, [T5J \13\ HI] ) .which 
first appeared in the context of black holes [TT] and later extended by Susskind [H] 
to string theory. According to the holographic principle, the entropy of a system 
scales not with its volume, but with its surface area. In the cosmological context, the 
holographic principle will set an upper bound on the entropy of the universe [15] . In 
the work [13], it was suggested that in quantum field theory a short distance cut-off 
is related to a long distance cut-off (infra-red cut-off L) due to the limit set by black 
hole formation, namely, if is the quantum zero-point energy density caused by a short 
distance cut-off, the total energy in a region of size L should not exceed the mass of 
a black hole of the same size, thus L 3 p^ < LM 2 . Applied to the dark energy issue, 
if we take the whole universe into account, then the vacuum energy related to this 
holographic principle is viewed as dark energy, usually called holographic dark en- 
ergy [13] [16], [17]. The largest L allowed is the one saturating this inequality so that 
we get the holographic dark energy density pa = 3c 2 M 2 L -2 where c 2 is a numerical 
constant and M~ 2 = 8nG. 

Choosing the Hubble horizon H~ l as the infrared cut-off, the resulting p\ is com- 
parable to the observational density of dark energy [IS], [IS]- However, in [I6J it was 
pointed out that in this case the resulting equation-of state parameter (EoS) is equal 
to zero, behaving as pressureless matter which cannot give accelerated expansion. 
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The particle horizon [T7] also results with an EoS parameter larger than —1/3, which 
is not enough to satisfy the current observational data, but the infrared cut-off given 
by the future event horizon [17], yields the desired result of accelerated expansion 
with an EoS parameter less than —1/3, despite the fact that it has problems with the 
causality. Another holographic DE model have been considered in [19], [20] . 
Based on dimensional arguments and the fact that the time derivative of the Hubble 
parameter naturally appears in the cosmological equations, in [21] we have proposed 
an infrared cut-off for the holographic density of the form p « aH 2 + (3H . Though the 
theoretical root of the holographic dark energy is still unknown, this proposal, which 
depends on the scale factor and its derivatives, may point in the correct direction as 
it can describe the dynamics of the late time cosmological evolution in a good agree- 
ment with the astrophysical observations. Another interesting fact of this model is 
that the resulting Hubble parameter (and hence the total density) contains a matter 
and radiation component [2T], which become relevant at high redshifts (radiation) 
in good agreement with the BBN theory, and (the matter component) explains the 
cosmic coincidence. An important fact is that this model can exhibit quintom nature 
without the need to introduce any exotic matter. This model has also proven to be 
useful in the reconstruction of different scalar field models of dark energy which repro- 
duce the late time cosmological dynamics in a way consistent with the observations 

[2211231 EI| ■ 

In this paper we study the cosmological constraints on the holographic dark energy 
model given in [2T] using the new dataset of type la supernovae (SNIa), called Union 
compilation. The Union compilation contains 414 SNIa and reduces to 307 SNIa after 
selection cuts. This 307 SNIa data set complemented with the CMB anisotropy and 
BAO (baryon acoustic oscillation) data, will be used to fit the parameter (3 of the 
model. Once (3 is known, it in turn fixes the parameter a and the integration constant 
via the flatness condition and the current (x=0) equation of state for the dark energy 



(see Eqs. (2.8) and (2.9) bellow). In the section 2 we review the main aspects of the 
model, in section 3 we fit the (3 parameter with the 307 SNIa data set, in section 4 
we use the joint SNIa, CMB and BAO data analysis to constraint (3 and show the 
best-fit (3 for each studied data set. The montecarlo simulation is used in section 
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5 to constraint the three parameters Q, m o,uo,f3 using the joint SNIa+CMB+BAO 
analysis, and in section 6 we give the statefinder and Om diagnosis to contrast the 
studied holographic DE model with the ACDM and other dynamical DE models. 



2 The Model 

Let us start with the main features of the holographic dark energy model. The 
holographic dark energy density is given by 

PA = 3(aH 2 + pH^j (2.1) 

where a and (3 are constants to be determined and H = a/a is the Hubble parameter. 
The usual Friedmann equation is 

H 2 = ^(p m + Pr + p A ) (2.2) 

where we have taken 8ttG = 1 and p m , p r terms are the contributions of non- 
relativistic matter and radiation, respectively. Setting x — In a, we can rewrite the 
Friedmann equation as follows 

H 2 = I (p m0 e- 3 * + Pr0 e-^) +*H 2 + ^ (2.3) 

Introducing the scaled Hubble expansion rate H = H/Hq, where Hq is the present 
value of the Hubble parameter (for x = 0), the above Friedman equation becomes 

H 2 = n m0 e- 3x + n r0 e- A * + aH 2 + ^ d ^- (2.4) 
where Q m o = p m o/?>Hl and Q r o = p r o/?>Hl are the current density parameters of 



non-relativistic matter and radiation. Solving the equation (2.4), we obtain 



where C is an integration constant. From this equation, the following equation for 
the holographic density is obtained 
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with pa = 
pressure 



^2 • The energy conservation equation gives the corresponding holographic 



PA 



2a - 3/3 - 2 



Ce -2x(a-l)/p + 



2(3 -a 



-Us 



(2.7) 



3/3 3(a- 2/3-1) 

There are three constants a, (3 and C which are related by the following two condi- 
tions. The first one is the restriction imposed by the flatness condition and can be 



obtained from Eq. (2.5) at x = 0: 

2 



-tt r0 + C = 1 



(2.8) 



3/3 -2a + 2 " tu 2/3-a + l 
The second condition is obtained by considering the equation of state for the present 
epoch values of the density and pressure (i.e. at x=0) of the dark energy p\ = cjo^ao 
(note that p\ = fl\o) 



2a - 3/3 - 2 

3/3 



C 



2(3 -a 



3(a- 2/3-1) 



n 



rO 



3(3 -2a 
2a -3/3- 2 



-n 



mO 



2/3 -a 
a -2(3-1 



f2 r n + C 



(2.9) 



Solving Eqs. (|2.8|) and (2.9) with respect to (3 it is obtained 

[2(1 - n m0 - n r0 ) + (3(n r0 + 3^(1 - n m0 - n r0 ) + 3)] (2.10) 



1 



a 



and 



C=l- 



2Q 



m 



2(^m0 + ^r()) — /3 [^rO + 3o;q(1 — f2 m — ^ro)] 

2^2 r Q 

2(fim0 + ^ro) — /5 [^rO + 3a>o(l — fi m Q — ^rd) — 1] 



(2.11) 



Replacing this expressions for a and C in (2.5) we obtain a /3-dependent Hubble 
parameter and cosmological evolution for a given matter density parameter Q mQ and 
present dark energy EOs parameter uj q . In previous works ( [221 El] ) we illustrated 
the behavior of the model by choosing some representative values for a and (3 guided 
by the redshift transition in the deceleration parameter and found the evolution of 
equation of state u(z). In general the obtained cosmological dynamics was very close 
to what is expected from the current astrophysical observations. We turn now to use 
several data sets to constrain the parameter (3 (and therefore a) of the holographic 



dark energy model (2.1). The constraints on the parameters of the holographic Ricci 
dark energy have been performed in 
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3 Constrains from SNIa observations 



We used the 307 super nova SNIa data from the union compilation set (table 11 from 
[2]). This 307 data are obtained after the cuts imposed on the total of 414 SNIa data 
using criteria of quality (2]. The data set gives the distance modulus fiobs(zi), and the 
theoretical (for a given model) distance modulus is defined by 



Vth(zi) = 5Log 10 D L (zi) + fi 



(3.1) 



where fio = 42.38 — 5Logioh, h is the Hubble constant Hq in units of 100 km/s/Mpc 
and Dl(z) = H cIl{z)/c. The luminosity distance times H is given by 

" 2 cdz' 



d L (z) = (l + z) 



o H(z' 



(3.2) 



where H(z, 9) from Eq. (2.5) in terms of z is given by 
2 



Hz, 



3(3 -2a + 2 



Q mQ (l + zY + 



1 



2(3 - a + 1 



O r0 (l + z) 4 + C(l + z) 



2(q~1)//3 



1 V2 



(3.3) 

with 6 = ((3, Q m , lu ) (after replacing a and C from ( 2.10|2.11 ) with f2 r0 = 0), but we 
will choose the values of fl m and Uq and constraint the constant (3. Next we minimize 
the statistical x 2 function (which determines likelihood function of the parameters) 
of the model parameters for the SNIa data. 



^2 ^ __ (Pobs (Zi) — Hth(Zi)) 2 



8=1 



(3.4) 



The x 2 function can be minimized with respect to the fiQ parameter, as it is indepen- 



dent of the data points and the data set [26]. Expanding the Eq. (3.4) with respect 
to no yields 

X 2 sn (6) = A(6)-2hoB(6)+h 2 C (3.5) 
which has a minimum for fio — B(9)/C, giving 

B{9f 



X 2 SN,rmn(0) = XsnW) = A(6) ~ 



c 



(3.6) 
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with 

A(Q} — CMo6g_Cg») ~ ^th{Zj,^Q = 0)) 

i=l a * 



307 

_ ^ Hobs(Zi) - fXth(Zj,flQ = 0) ^ ^ 

i=l a i 



307 

= Ei 



Now this x|tv i s independent of /io and can be minimized with respect to the parame- 
ters of the theoretical model. In our analysis we will take some samples of Q m and co>o 
from the observational data and minimize with respect to the remaining parameter 
(3. The best fit value for (3 with la uncertainty and the corresponding x m im from the 
analysis of the Union sample of 307 SNIa pj is resumed in table I. 



^m0 


UJ 


(3{la) 


a 


h 


y 2 . 
Amm 


0.28 


-1.29 


292 +0 - 040 
u.zyz_ 035 


0.751 


0.708 


312.131 


0.28 


-1 


u.oot_ 106 


0.965 


0.698 


312.744 


0.28 


-0.969 


u.uoo_ 129 


1.017 


0.696 


313.189 


0.29 


-1.29 


n oi O+0.044 
U ' Jiz -0.038 


0.749 


0.707 


312.035 


0.29 


-1 


U.OOO_o.i20 


0.986 


0.697 


312.931 


0.29 


-0.969 


7+0 .214 
O-' 1 ' -0.146 


1.046 


0.696 


313.414 


0.3 


-1.29 


D QQQ+0.048 
u - ooo -0.041 


0.748 


0.707 


311.951 


0.3 


-1 


0.693^i| 


1.012 


0.697 


313.138 


0.3 


-0.969 


0.7871™ 


1.080 


0.695 


313.659 



Table 1: The best-fit values for (3 with la error, from the SN la analysis. 



The best fit value for h can be obtained from the relation fiQ = B((3)/C at the 
best fit value for (3. We used three representative values of Q m o combined with three 
representative values of Uq, obtaining a total of 9 best fit values for (3. Note that (3 
significantly changes with the change in uu and the best fit h is less sensitive, taking 
values in the region h ~ 0.7. Fig. 1 shows the behavior of the Hubble parameter, 
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using the SN la 307 data set, for three different best-fit values of (3 taken from table 
I, and the corresponding likelihood behavior for the same three values. 




Figure 1: The observational H(z) data with error bars, and the holographic model 
H(z) for three different best fit values of (3 from the 307 SN la analysis, and (right) 
the corresponding likelihood behavior for the same set of parameters. 

Note from table I, that in all cases for ujq > —1, the constant a > 1 giving a 
quintessence character to the model, as the power of (1 + z) in the last term in Eq. 



3.3 becomes positive; and in the cases with lj < —1, a changes to a < 1 and the 



power of (1 + z) becomes negative giving a quintom-like character to the model. 
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4 Constraining from the SNIa+CMB+BAO 



In this section we constrain (3 in the holographic dark energy model (2.1) by using 
the SNIa combined with the CMB and BAO observational data. The shift parameter 
R [27] from the cosmic microwave background (CMB) anisotropy, and the distance 
parameter A of the measurement of the baryon acoustic oscillation (BAO) peak in the 
distribution of SDSS luminous red galaxies |4j, are also used extensively in obtaining 
the cosmological constraints. The shift parameter R is defined by 

R = / S- (4.1) 
m °J H(z) 1 ' 

where z = 1090 is the redshift of the recombination [28]. The distance parameter A 
is given by [2H] 

"i r» dz 

_z b Jo H{z) 

with zt = 0.35. And now, using the combined data of the 307 Union SN la, the shift 
parameter R of CMB and the distance parameter A of BAO, we perform the joint 
analysis to constraint the constant (3. The total \ 2 is given by 

X 2 = X 2 SN + X 2 cmmb + X 2 bao ( 4 - 3 ) 
The best-fit model parameter can be determined by minimizing the total \ 2 ■ Here 

X 2 CMB and X 2 BAO are S iven b y 

2 _ (R — Rpbs) 2 2 _ (A ~ Apbs) 2 ,x 

XCMB — 2 ' XBAO — 2 K^-^) 



A = rt'lHizb)- 11 * 



2/3 



(4.2) 



The SDSS BAO measurement ([2H]) gives the observed value of 
A = 0.469(n s /0.98)" 35 ± 0.017 with the spectral index n s as measured by WMAP5 
[28], taken to be n s = 0.960. The value of the shift parameter R has also been up- 
dated by WMAP5 [2H] to be 1.710 ±0.019. The table II shows the best fit value for (3 



with la uncertainty and the corresponding x m ini fr° m the joint analysis of the Union 
sample of 307 SNIa, CMB and BAO observations. We used a combination of three 
different values for Q m0 and uq respectively. 
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f2 m 


UJ 


(3{la) 


a 


h 


y 2 . 

A.min 


0.28 


-1.29 


0.437 + °;°l° 


0.767 


0.718 


321.028 


0.28 


-1 


0.625_q 


0.983 


0.699 


312.828 


0.28 


-0.969 


0.6541^24 


1.017 


0.696 


313.196 


0.29 


-1.29 


0.433!{ffiS 


0.765 


0.715 


319.404 


0.29 


-1 


0.618lg£t 


0.979 


0.697 


313.134 


0.29 


-0.969 


0.646j£™ 


1.012 


0.695 


313.663 


0.3 


-1.29 


n 43n+o-023 

u -^ ou -0.016 


0.763 


0.713 


318.597 


0.3 


-1 




0.975 


0.695 


314.055 


0.3 


-0.969 


0.638« 


1.008 


0.693 


314.727 



Table 2: The best-fit values for (3 with la error, from the joint SNIa+CMB+BAO 
analysis. 



Note that R and A are independent of /j,q, but the best fit value for h is obtained 
from fiQ = B(/3)/C at the new best-fit values for (3. Fig. 2 shows the behavior of H(z) 
for the joint analysis of the SNIa+CMB+BAO data, and the likelihood behavior for 
the SNIa+CMB+BAO analysis, for three representative best-fit values of f3 taken 
from table II. 

Looking at tables I and II, and analyzing all the obtained values for x m im the best 
fit for the SN la analysis happens at (3 — 0.333^q;q4^, corresponding to the lowest 
= 311.951 with Q m0 = 0.3 and Uq = 



analysis, at (3 — 0.625^n°J3, lowest xL 



1.29, and for the joint SN Ia+CMB+BAO 
0.023' """a™ Amin ~= 312.828 with tt m0 = 0.28 and lu = -1. 
Fig. 3 shows the behavior of the Hubble parameter with the redshift for the lowest 
X m in from the SN la analysis and the SNIa+CMB+BAO joint analysis, and the 
corresponding likelihood behavior. 

Note that we have constrained (3 alone, fixing each time the current values of f2 m0 
and uiq, and hence we can not make definite conclusion about the quintessence or 



quintom nature of the model (2.1). Nevertheless looking at Fig. 3 for experimental 
(error bars) versus holographic model H(z), it can be seen that the quintom behavior 
(a < 1) is favored by both, the SNIa and joint SNIa+CMB+BAO analysis. 
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Figure 2: The observational H(z) data with error bars, and the holographic model 
H(z) for three different best fit values of (3 from the joint 307 SNIa+CMB+BAO 
analysis, and (right) the corresponding likelihood behavior for the same set of param- 
eters. 




0.0 0.5 1.0 1.5 0.0 0.2 0.4 0.6 0.8 1.0 

z fS 

Figure 3: The observational H(z) data with error bars, and the holographic model 
H(z) for the best-fit f3 corresponding to the lower Xmin from the SNIa (table I) and the 
joint SNIa+CMB+BAO (table II) analysis, and (right) the corresponding likelihood 
behavior. 
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To better illustrate the cosmological dynamics, in Fig. 4 we plot the effective equation 
of state parameter (u e ff = Pa/(pa + p m ) and the deceleration parameter q(z) for the 
two best-fit (3 used in Fig. 3. 







SNIa + CMB + BAO 




- /J=0.625,n mO =0.28,« =-l 






/3=0.333, n m0 =0.3, wo=-1.29 - 




SNIa 



SNIa 

0.5 - ^=0.333, n m0 =0.3, in =-1.29, z r =0.45 




SNIa + CMB + BAO 
/3=0.625, n m0 =0.28, u„=-l, z r =0.69 



-0.5 0.0 



Figure 4: The the evolution of the effective equation- of -state u e ff(z), for the best-fit 
j3 from the SNIa (table I) and the joint SNIa+CMB+BAO (table II) analysis, and 
(right) the deceleration parameter showing the redshift transition z? for both cases. 
The redshift transitions Zt are in the range accepted by observations 



5 Constraining the cosmological and model param- 
eters using Montecarlo 

In this section we perform constraints on the three parameters (3, Q m o and uo of our 
holographic dark energy model, combining the observations from 307 SNIa, CMB 
and BAO, by using the montecarlo technique to restrict the Xmin an d choosing an 
appropriate intervals for (3, Q m o and ujq. The analysis was done over the total of 10 6 
data, with the restrictions on Xmin an d P taken according to the previous fits consigned 
in tables I and II: xLn < 313.000, 0.25 < (3 < 2.3, 0.25 < Q m0 < 0.35 and -1.3 < 
ujq < —0.8. After 10 6 iterations, the best fit-value for the triplet of model parameters 
was (3 = 0.593^;^(1(7), tt m0 = 0.283ljffi(l<7) and lu = -1.036l^I(la), with 
Xmin = 312.734. The value of c^o is slightly lower than —1 showing that the model 
tends to behave as quintom. Fig. 5 shows the Hubble parameter versus z for the 
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best-fit values of (j3, Q m o, ^o), and the evolution of the effective and holographic EoS 
parameters. 



fl,„=0.283, £=0.593, w =- 1.036 




£=0.593, n m0 =0.283, u =-1.036 






"A 




SNIa + CMB + BAO, Monte Carlo - 




(j efec 



Figure 5: The observational H(z) data with error bars, and the holographic H(z) 
for the best-fit model parameters (f3,Q m0 ,uj ) . The right graphic shows the effective 
EoS and the holographic EoS u\(z). The transition deceleration acceleration (oo e ff — 
— 1/3,) occurs at zt ~ 0.67. The quintom character of the model with the u>\ = —1 
crossing is seen from the u\ behavior. 



6 The Statefinder and Om diagnostics 



To better understand the properties of the DE model (2.1) is useful to compare 
with the model independent diagnostics which are able to differentiate between a 
wide variety of dynamical DE models, including the ACDM model [5U]. We will use 
the diagnostic introduced in [30] known as statefinder, which introduces a pair of 
parameters (r, s) defined as 

r= ^' s= ^h) (61) 

which are proven to be useful to characterize a given DE model, as this pair depend 
only on high derivatives of the scale factor a and 'a, acquiring some geometrical sense 
as depend only on the spacetime metric. In terms of the total density and pressure 
(r, s) can be written as ([!]). 

r = l+ 9(p + P) g, s= iP±Ai (6 . 2) 
2p p p p 
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Is clearly seen from Eq. (6.2) that in the flat FRW background, the ACDM model 



corresponds to a fixed point (s = 0, r — 1) in the r — s plane. The trajectories in the 
s — r plane corresponding to different DE models may exhibit qualitatively different 
behaviors, and this is a good way to establish the departure of a given DE model 
from the ACDM. In this work we apply the statefinder diagnostic to the holographic 
DE model ( 2.1[ ) (|2"Tj). using the calculated above best-fit parameters with the SNIa 
and joint SNIa+CMB+BAO data analysis. 



The statefinder parameters for the model (2.1) are given by 

C(a- 1){2 - 2a + 3(3) 2 e 3x 



1 + 



and 



f3 2 [C(2a -3(3- 2)e 3x - 2Q m0 e 2x ( a -^ 



a 







(6.3) 



(6.4) 



Note that a and C have to be replaced by solutions (2.10) and (2.11). From this 
Equations one can see that s = and r = 1 for a = 1, which corresponds to the 
ACDM model. To illustrate our case, in Fig. 6 we plot the statefinder diagram in the 
s—r plane for the best-fit (/?, £l m o, ^o) taken from the joint SNIa+CMB+BAO analysis 
using montacarlo and for the best-fit (3 = 0.654, a = 1.017 (Q m o = 0.28, ojq = —0.969) 
taken from table II, for comparison. 



Today 



-0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 



Figure 6: The statefinder evolution for the holographic model in the s — r plane with 
((3 = 0.593, a = 0.946, C = 0.7, fi m0 = 0.283) (arrow up) and ((3 = 0.654, a = 1.017, 
C = 0.71, Q m0 = 0.28) (arrow down). 
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Since s is constant, the trajectory in the s — r plane is a vertical line with r 



monotonically increasing from r = 1 to r ~ 1.29 (at x— > oo in (6.3)) for the used 
best-fit values. The current (s, r) values are (s = — 0.06,r(x = 0) = 1.203). The 
SCDM (standard cold dark matter) model is shown in the point (s = l,r = 1). The 
evolution trajectory in the s — r plane is similar to the one obtained for the Ricci DE 
model [31] • Note that for a > 1 the r trajectory points in the opposite direction. 
Another useful statefinder diagram is given by the trajectory in the q — r plane. The 
deceleration parameter q 



1 

q= 2 

and r in terms of q is 



1 + 



applied to the model (|2.l|) is given by 

C(2 - 2a + 3(3) 2 e 3x 



(5 [C(2a -3(3- 2)e 3a; - 2Q m0 e 2 < a - 1 )IP]_ 



a 







(2q 



(6.5) 



(6.6) 



Fig. 7 shows the evolution trajectory of the model in the q — r plane for the same 
best-fit values of fig. 6. 



Today 



-1.0 -O.S -0.6 -0.4 -0.2 0.0 0.2 0.4 



Figure 7: The statefinder evolution for the holographic model in the q — r plane, for 
(p = 0.593, a = 0.946, C = 0.7, tt m0 = 0.283) and (/? = 0.654, a = 1.017, C = 0.71, 
n m0 = 0.28). 

we see that the ACDM and the holographic model start diverging from the same 
point in the past (the standard cold dark matter SCDM). The current values of the 
statefinder parameters are in q = —0.614 and r = 1.203. This trajectory has constant 
negative slope, corresponding to a quintom behavior of the model with a < 1. The 
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trajectory for a = 1.07 taken from table II (for illustration) have positive slope as 



can also be seen from Eq. (6.6) for a > 1, characterizing the quintessence behavior 



of the proposed holographic DE model. In this plane the ACDM scenario is the 



horizontal line in Fig. 7 (as follows from Eq. (6.6)) , and evolves from the SCDM 
in the past (q = 0.5, r = 1), ending at the steady state cosmology (SS) in the future 
(q = -l, r = 1). 

An alternative way of distinguishing ACDM from other DE models without directly 
involving the EoS, is the Om diagnostic introduced in [32]. This Om diagnostic is 
constructed from the Hubble parameter, depending only on the first derivative of 
the expansion factor a(t), and hence depends only upon the expansion history of our 
Universe. It is defined by 

H 2 (x) - 1 

Om(x) = _\ X J - (6.7) 



1 



which for the present holographic model turns out to be 

Om(x) 



2Qmo r -3x I ry p -2x(a-l)/f3 _ i 

2-2q+3/3 c -T x 



-3x 



(6i 



for DE models with constant EoS the Om(x) takes a simple form (see [32]) and for 
the ACDM model takes the simplest form, being equal to the density parameter [32] . 
In Fig. 8 we show the Om diagnostic for the holographic model with the best-fit 
(/?, Qrno, ujq) taken from the joint SNIa+CMB+BAO analysis using montacarlo, and 
for the best-fit (3 = 0.654, a = 1.017 {Q m0 = 0.28, u = -0.969) taken from table II 
for comparison. 
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Figure 8: The Om(x) for the holographic model with the montecarlo best-fit (3 = 0.593, 
a = 0.946, C = 0.7 and Q m0 = 0.283, and for the best-fit (3 = 0.654, a = 1.017 from 
table II. 

In this graphic the evolutionary difference between the ACDM and the studied 
holographic DE model is more clear, with the Om curve turning down for a < 1 and 
up for a > 1. 



7 discussion 

In this letter we have obtained the constraints on the parameters of the holographic 



dark energy model described by the density Eq. (2.1) [21], using the latest observa- 
tional data including the 307 Union sample of SNIa, the CMB shift parameter given 
by WMAP5, and BAO measurement from SDSS. We first used the 307 Union sample 
of SNIa to constraint the parameter (3, assuming priors about fundamental cosmo- 
logical quantities such as the current EoS for dark energy and the matter density 
parameter, and have generated the table I with different priors proposed. Looking 
at the minimum of the different Xmin listed in table I, a representative best-fit can 
be given by (3 = 0.333+o;q4i (la) with priors Vt m o = 0.3 and ujq = —1.29. The same 
considerations for the joint SNIa+CMB+BAO analysis resumed in table II, give the 
best-fit (3 = 0.625^qq23 with priors fl m o = 0.28 and u = —1, with la uncertainty. In 
fig. 3 we plot the evolution of H(z) for the two best-fit j3 for each analysis. Of course, 
the assumed priors may subvert the efficacy of the method and would be better to 
minimize the x 2 with respect to all the relevant cosmological parameters. In this 
sense we considered useful to use the montecarlo method to generate the Xmin with- 
out assuming priors about the relevant parameters, but restricting their range to a 
reasonable intervals (between the limits on Q m o and lvq set by experiments). After 10 6 
iterations we obtained the following best fits with la uncertainty: (3 = 0.593^q;q23, 
tt m0 = 0.283+°;°};j and cu = -1.036+°;°^ ; with X 2 min = 312.734. The value of u 
states that the possibility of quintom behavior can not be discounted. The behavior 
of H(z) contrasted with the observational error bars is shown in Fig. 5. Note also 
that the values of best-fit h were found in a very narrow intervals (see tables I and 
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II): for the 307 SNIa data 0.695 < h < 0.708 and for the joint SNIa+CMB+BAO 
0.697 < h < 0.718, which are in the limits set by different experiments [3J, [SJ, [33J. 
From Figs. 3 and 5 for H(z) and using the standard deviation from the observational 
error bars, the best-fit curve using the standard deviation criteria was obtained for 
the montacarlo SNIa+CMB+BAO analysis (Fig. 5). 

we also used the statefinder and Om diagnosis to differentiate this model with other 
models of dark energy. With the statefinder diagnosis the evolutionary trajectory 
in both, the s — r and q — r planes points in opposite directions depending on the 
quintessence (a > 1) or quintom (a < 1) nature of the model. In the q — r plane the 
trajectory diverges from the SCDM point, but never reaches the ACDM line (except 
for a = 1), making the difference with other DE models [HD], [34]. With the Om 
diagnosis there is a marked difference between the holographic model and the ACDM 
at the current time, and in the case of a < 1 the holographic trajectory intercepts 
in the past (at z ~ 0.65) the ACDM line. Finally we hope that future high precision 
experiments allow accurately determine the parameters of the model to define its 
quintessence or quintom nature. 
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